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Abstract — A global model is presented that can be used to 
study attitude maneuvers of a rigid spacecraft in a circular 
orbit about a large central body. The model includes gravity 
gradient effects that arise from the non-uniform gravity field 
and characterizes the spacecraft attitude with respect to the 
uniformly rotating local vertical local horizontal coordinate 
frame. An accurate computational approach for solving a 
nonlinear boundary value problem is proposed, assuming that 
control torque impulses can be applied at initiation and at 
termination of the maneuver. If the terminal attitude condition 
is relaxed, then an accurate computational approach for solving 
the minimal impulse optimal control problem is presented. 
Since the attitude is represented by a rotation matrix, this ap- 
proach avoids any singularity or ambiguity arising from other 
attitude representations such as Euler angles or quaternions. 

I. Introduction 

The attitude dynamics of an uncontrolled rigid spacecraft 
in a circular orbit about a large central body, including 
gravity gradient effects, have been extensively studied; see 
[1], [2]. There are 24 distinct relative equilibria for which the 
principal axes are exactly aligned with the local vertical local 
horizontal (LVLH) axes, and the spacecraft angular velocity 
is identical to the orbital angular velocity of the LVLH 
coordinate frame. Linear rotational equations of motion that 
describe small perturbations from any relative equilibrium 
solutions are well known. Linear attitude control of a rigid 
spacecraft in a circular orbit, including linear gravity gradient 
effects, has also been addressed in [2]. However, linear 
controllers have the limitation that they are only applicable 
to small attitude change maneuvers. 

The emphasis in this paper is on large angle attitude ma- 
neuvers of a rigid spacecraft in a circular orbit about a large 
central body, including gravity gradient effects. A nonlinear, 
globally defined model is introduced. This model describes 
the attitude of the spacecraft, relative to the uniformly 
rotating LVLH coordinate frame, by a rotation matrix. The 
model includes gravity gradient terms that reflect the rotation 
of the LVLH frame, and terms that reflect control input 
torques. In the problems studied in this paper, independent 
impulsive control torques can be applied about each principal 
axis. Gravity gradient moments are significant in Earth orbits 
for orbit altitudes between 400 km and 40,000 km and for 
attitude maneuver times that are not small compared with 
the orbital period. 
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LVLH frame 

Fig. 1. Coordinate frames 

We study open loop attitude maneuvers that can be ac- 
complished by using two impulsive torque controls, one 
occurring at the initial time and one occurring at the terminal 
time of the maneuver. The attitude motion in between the 
initial time and the final time is uncontrolled. Two classes 
of attitude maneuver problems are studied. In Section [ill] a 
terminal attitude is specified so that there is a unique impulse 
sequence satisfying the boundary conditions. This problem 
can be solved computationally use a root finding algorithm. 
In Section IIVI the specified terminal attitude condition is 
relaxed such that there are many impulse sequences that 
satisfy the boundary conditions. In this case we seek the 
minimum total impulse sequence that satisfies the boundary 
conditions. This problem can be solved computationally 
using a constrained minimization algorithm. 

For both classes of attitude maneuver problems, our 
contribution is to demonstrate how sensitivity derivatives, 
used in the computational algorithms, can be determined 
effectively and accurately such that they satisfy the global 
geometry of the problem. Since the attitude is represented by 
a rotation matrix in the special orthogonal group SO (3), and 
the sensitivity derivatives are expressed in terms of the Lie 
algebra so(3), this approach completely avoids singularities 
or ambiguities that arise from other representations of the 
rotation group, such as Euler angles and quaternions. 

II. Rigid spacecraft models in a circular orbit 

We assume that a rigid spacecraft is on a circular orbit with 
a constant orbital angular velocity uio e K. In this section, the 
continuous equations of motion and a geometric numerical 
integrator, referred to as Lie group variational integrator, for 
the attitude maneuver of a spacecraft in a circular orbit are 
given following the results of [3], We identify TSO(3) ~ 



S0(3) xso(3) by left translation, and we identify so (3) ~ M. 3 
by an isomorphism S(-) : R 3 i— > so(3). 

Continuous equations of motion: We define three rotation 
matrices in SO(3); 

R : from the body fixed frame to the inertial frame, 
R : from the LVLH frame to the inertial frame, 
R bl : from the body fixed frame to the LVLH frame, 

where the inertial frame and the LVLH frame are illustrated 
in Fig.Q] Thus, R bl = R uT R w . 

The on-orbit spacecraft equations of motion are given by 
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where n, il G R 3 are the angular momentum and the angular 
velocity of the spacecraft expressed in the body fixed frame, 
respectively, and M 9 G R 3 is the gravity gradient moment. 
The isomorphism between R 3 and so(3) is defined such 
that S{x)y = x x y for any x, y G R 3 . Since the orbital 
angular velocity w is constant, the solution of (|3} is given 
by R u (t) = R H (0)e s ^ oe ^ t . 

Gravity gradient moment: The gravity gradient moment 
is derived in [2]. We present an alternative way to obtain the 
gravity gradient moment directly using the gravity potential; 

GM 
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where x G R 3 is the position of the spacecraft in the inertial 
frame, and p G R 3 is a vector from the center of mass of the 
spacecraft to a mass element in the body fixed frame. G is 
the gravitational constant and M is the mass of the Earth. 

^From [3], the gravity gradient moment M 9 can be 
determined by using the following relationship; 
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We derive a closed form for M 9 from (0, by assuming that 
the spacecraft is on a circular orbit so that the norm of x is 
constant, and the size of the spacecraft is much smaller than 
the size of the orbit. 

As shown in Fig. [2 the coordinate of the spacecraft 
position in the LVLH frame is roe 3 , where ro G R is the 
radius of the circular orbit, and e 3 = [0, 0, 1] T . Therefore, 
the position of the spacecraft in the inertial frame is given 
by x — roR h e 3 . Using this expression, 
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where p = tAt G R 3 is the unit vector along the direction 
of p. Assume that the size of the spacecraft is significantly 



smaller than the size of the orbit, i.e. ^ <C 1. Using a Tay- 
lor series expansion, we obtain the 2nd order approximation. 
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Since the body fixed frame is located at the mass center 
of the spacecraft, j B p dm = 0. Therefore, the first term 
in the above equation vanishes. Because e^R bl p is a scalar 
quantity, we can rewrite the above equation as 
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where ujo = \/ G R is the orbital angular velocity, and 

V r o 

J G R 3x3 is the moment of inertia matrix of the spacecraft. 
Substituting (|6jl into (|5jl, and using the property S(x x y) = 
yx T — xy T for x, y G R 3 , we obtain an expression for the 
gravity gradient moment as follows. 

M 9 = 3ulR blT e 3 x JR blT e 3 . (7) 

Discrete equations of motion: In the continuous equations 
of motion, the structure of (|2ji, 0, and @ ensures that 
Rr, R h , and R bl evolve on the special orthogonal group, 
SO (3). However, general numerical integration methods, 
including the popular Runge-Kutta methods, do not preserve 
the orthogonality property of this group. For example, if we 
integrate @ by a typical Runge-Kutta scheme, the quantity 
R bl R bl inevitably drifts from the identity matrix as the 
simulation time increases. 

The rotation matrix is commonly parameterized by Euler 
angles or quaternions. These attitude kinematics equations 
can be numerically integrated and are used to recompute 
the rotation matrix. However, Euler angles are not global 
expressions of the attitude since they have associated singu- 
larities. The analytical expressions for sensitivities are hard 
to develop since many trigonometric terms are encountered. 
Quaternions have no singularity, but quaternions must lie on 
the three sphere S 3 . General numerical integration methods 
do not preserve the unit length of a quaternion. There- 
fore, quaternions have the same numerical drift problem 
as rotation matrices. Furthermore, quaternions, which are 
diffeomorphic to SU(2), double covers SO(3). So there are 
inevitable ambiguities in expressing the attitude. 

These cause significant inaccuracies in numerical simula- 
tions based on quaternions and Euler angles. In particular, 
the gravity gradient moment, as given in @, depends on 
R u directly, and consequently, errors in computing R bl cause 
errors in the gravity gradient moment. These effects are more 
pronounced when the simulation time is large. 

Lie group variational integrators preserve the orthonormal 
structure of SO (3) without any reprojection or parametriza- 
tion. They also conserve the momentum map, and the sym- 
plectic property of rigid body dynamics. In addition, the total 
energy is well-behaved, as it only oscillates in a bounded 
fashion about its true value. So, Lie group variational inte- 
grators are geometrically exact. Using the results given in [3], 



a Lie group variational integrator for the attitude dynamics 
of a spacecraft in a circular orbit are given by 
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where the subscript k denotes variables at the fcth time step, 
and h € R is the integration step size. The matrix Jd £ R 3x3 
is a nonstandard moment of inertia matrix defined by Jd = 
±tr[J]/ 3x3 - J. 

The matrix F k = R b f R b k l +l e SO (3) is the relative 
attitude between integration steps, and it is obtained by 
solving (gj. Since F k and e s ^° e ^ h are in SO(3), and SO(3) 
is closed under matrix multiplication, R°}, R l k , and R b k 
evolve in SO (3) automatically for all k according to dlOi . 
(lilt , and (I12> . The actual computation of F k is done in the 
Lie algebra so (3) of dimension 3, and the rotation matrices 
are updated by multiplication with the exponential of a skew- 
symmetric matrix. So, Lie group variational integrators are 
numerically efficient, and there is no excessive computational 
burden in updating the 9 elements of the rotation matrix. The 
properties of these discrete equations of motion are discussed 
more explicitly in [3] and [4]. 

Since ujq is constant, the solution of il Q is given by 



III. Spacecraft attitude maneuvers 
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A. Problem formulation 

A two point boundary value problem is formulated for 
spacecraft rest-to-rest maneuver between two given attitudes 
for a fixed maneuver time. This is a Lambert type boundary 
value problem on SO(3). The initial attitude and the desired 
terminal attitude are expressed as rotation matrices with 
respect to the LVLH frame, namely Rff, R^ e SO(3). Two 
impulsive control torques are applied at the initial time and 
the terminal time. We assume that the control torques are 
purely impulsive, which means that each impulse changes 
the angular velocity of the spacecraft instantaneously, but it 
does not have any effect on the attitude of the spacecraft at 
that instant. The rotational motion of the spacecraft between 
the initial time and the terminal time is uncontrolled. This 
assumption is reasonable when the operating time of the 
spacecraft control moment actuator is much smaller than the 
fixed maneuver time. 

We define this boundary value problem directly on SO (3) 
instead of using parameterizations such as Euler angles and 
quaternions. This approach allows us to define the attitude of 
the spacecraft globally without singularities and ambiguities. 
We use the discrete equations of motion given in d8l- J12> 
for the problem formulation and for the following analysis. 

We transform this two point boundary value problem into 
a nonlinear root finding problem. For a given initial angular 



momentum IIo and an initial attitude of the spacecraft Rq 1 , 
the terminal angular momentum II jy and the terminal attitude 
are determined by the discrete equations of motion. By 
choosing the initial angular momentum so that the terminal 
attitude of the spacecraft is equal to the desired attitude, i.e. 
Rjy = R b N d , we obtain the initial impulse and the terminal 
impulse. Thus, the nonlinear boundary value problem for the 
spacecraft attitude maneuver is formulated as 

given :R bl ,R b 'N 



such that R 
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where N S N is the number of integration steps determined 
by N — j- for the fixed maneuver time T and the fixed 
integration step size h. 

B. Computational approach 

We solve a sequence of linear boundary value problems 
whose solutions converge to the solution of the nonlinear 
boundary value problem. 

Linearization: The equations of motion are linearized 
about a given trajectory, and they are expressed in terms 
of the Lie algebra so (3). Consider small perturbations from 
a given trajectory denoted by II|, R k ' e ,R k ' e , F k . 
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where e 6 R. Since the orbital angular velocity is constant, 
8Rl =0. 

The infinitesimal variation of the angular momentum 5Ii k 
can be expressed in R 3 . The variation of the rotation matrix 
R" k %£ e SO (3) can be expressed as 

R%>* = R%e' s ^\ 

where ( k £ R 3 and S(( k ) £ so(3) is a skew-symmetric 
matrix. Since the map S(-) is an isomorphism between so(3) 
and M 3 , ( k is well defined. Then, the infinitesimal variation 
5R b k is given by 
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SF k is obtained from definition fllOi , and j!8l > as 
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In summary, equations dl8> . dl9l and l|2(J 
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Substituting (TBI . ( fTBT l. and jTTJ into (|8}, and and 

ignoring higher-order terms, the linearized discrete equations 
of motion are given by 

sn k+1 = SFfu k + F^sn k 

h „ „ h rp „ „ h „ „ (21) 
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These equations are not in standard form since 1221 is an 
implicit equation in Ji 7 ^. By using ( 1181 . ( I19> and ( I20L we 
will obtain an explicit solution of J22I . and rewrite the above 
equations in standard form. 

Substituting ( I20t into ( I23l l. and using the property 
S(x)y = —S(y)x for all x,y G R 3 , 6M^ can be written 
as 
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where .M fc G R 3x3 . 

Using ([19} and the property S(Rx) = RS(x)R T for x G 
M 3 , i? G SO (3), the right hand side of {23 is written as 

5F k J d - JdSF? 

= -{S(Ck)F k J d + J d F k T S(Ck)} 
+ {S(F k C k +i)F k J d + J d F£S(F k <; k+1 )} . (25) 

Substituting 1241 and d25i into 1221 . and using the property 

S{x)A + A T S(x) = S({ti[A] hxs -A}x) for x G M 3 , 
A G R 3x3 , (I22> can be transformed into an equivalent vector 
form; 

7i 2 
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Multiplying both sides by F^ {ti[F k J d ] I 3x3 - F fe J d } _1 and 
rearranging, we obtain 
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where .4*, G 
obtained as 



p3x3 



Substituting {26j into (JT9J, 5F k is 



d"-Ffc = -7?fe-Ffe + F k r)k+i, 

= -S((k)Fk + FkS{(k+i), 

= S((k)Fk + F k S{Ak(k + BkSUk) 

Equation i21l is an explicit solution of (122b . 
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Now we rewrite (12 U in the standard form of linear discrete 
equations of motion. Substituting (I24i . J27i into ( I21t . and 
rearranging, we obtain 
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where C k ,T>k G R 3x3 . 

In summary, (I26i and J28t are linear discrete equations 
equivalent to (12 It . J22i and ( I23l >. and they can be written as 
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where A k G R 6x6 . Equation d29t is the linear discrete 
equation for perturbations of the attitude dynamics of a 
spacecraft in a circular orbit, expressed in terms of K 3 ~ 
so (3). The important feature of d29t is that it is linearized 
in such a way that it respects the geometry of the special 
orthogonal group SO(3). 

Linear boundary value problem: The solution of \29\ is 
given by 
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where $^ G 

For the given boundary value problem, Co = since the 
initial attitude of the spacecraft is given and fixed, and SUn 
is free since the terminal angular momentum is compensated 
by the terminal impulse. Then, we obtain 
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This equation provides $12, the sensitivity derivative of 
the terminal attitude with respect to a change in the initial 
angular momentum. It states that for a given trajectory, if we 
update the initial angular velocity by JIIo, then the terminal 
attitude is changed from R% to i?^e 5 ^ N ) = R^e 8 ^ 12511 ^. 

We choose the change of the initial angular momentum 
so that the updated terminal attitude is equal to the desired 
terminal attitude; R'$e s( - iS>12SlUl ' > = R^ , or equivalently, 

dn a = $r 2 1 S'~ 1 (logm(< T < £i )) , (31) 

where S f_1 (-) : So(3) 1— > R 3 is the inverse mapping of 
S(-), and logm denotes the matrix logarithm. Equation ( 13 11 1 
provides a solution of the linear boundary value problem 
for the attitude dynamics of a spacecraft in a circular orbit, 
assuming that $12 is invertible. 



Nonlinear boundary value problem: The linear boundary 
value problem is solved successively so that its solution 
converges to the solution of the nonlinear boundary value 
problem. A numerical algorithm is summarized as follows. 



Set Error = 2e s . 
Guess an initial condition II 
Set i = 0. 
while Error > £5. 
Find njj*\i# W 
Compute the error; 
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Update the initial condition; II 
Set i = i + l. 
end while 
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Here the superscript (i) denotes the ith iteration, and es, c £ 
R are a stopping criterion and a scaling factor, respectively. 

This computational approach utilizes an exact and efficient 
method to compute sensitivity derivatives in the special 
orthogonal group SO(3). The sensitivity derivatives are then 
used to solve the two point boundary value problem. 

C. Numerical example 

Three spacecraft rest-to-rest maneuvers between relative 
equilibrium attitudes are considered. The resulting motions 
are highly nonlinear, large angle maneuvers. 

The mass, length and time dimensions are normalized. 
The moment of inertia of the spacecraft and simulation 
parameters are chosen as J = diag [1, 2.8, 2], e s = 10~ 14 , 
c = 0.1, h = 0.001. Each maneuver is completed in a 
quarter of the orbit, T = ^ . The boundary conditions and the 
corresponding computed impulsive control are as follows. 

(i) Rotational maneuver about the LVLH axis e\\ 



i?o=/ 3 x 3 , R"n = diag [1,-1,-1] : 



IE, = [2.116, 1.531, -1.782] J , 
II 7v = [-2.116, 1.531, 1.782] T . 
(ii) Rotational maneuver about the LVLH axes e\ and e 2 - 



R b l = diag [1,-1,-1] 



R bl 
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— n^y = [0.397, -1.586, -1.310^ . 
(iii) Rotational maneuver about the LVLH axes e 2 and 63: 



R b l = diag [1,-1,-1] 



nbl _ 
R N d ~ 





-1 


\T 



n = [1.047, 0.437, 2.800] 
n^y = [-1.416, -1.761, 1.159] T . 

Fig. [2] shows the attitude maneuver of the spacecraft, 
and the angular velocity response for each case. (Simple 
animations which show these spacecraft maneuvers can be 
found at |http : / / www . umich . edu/ ~tylee| ) 



IV. Optimal spacecraft attitude maneuvers 

A. Problem formulation 

An optimization problem is formulated as a rest-to-rest 
maneuver of an axially-symmetric spacecraft from a given 
initial attitude to a given terminal reduced attitude for a fixed 
maneuver time. The initial attitude is expressed by a rotation 
matrix with respect to the LVLH frame, namely Rfj. The 
terminal desired attitude is given by the reduced attitude, 
Ajv d = -Rjv e 3 £ § 2 . This reduced attitude represents the 
direction of the spacecraft axis of symmetry 63 in the LVLH 
frame. Two impulsive control moments are applied at the 
initial time and the terminal time, and the maneuver of the 
spacecraft between the initial time and the terminal time is 
uncontrolled. 

In the problem studied in section Qll] the initial parameter 
n is exactly prescribed by the constraint R!f] d = R%, and 
the discrete dynamics. In this section, we relax the terminal 
constraint by only specifying it up to a rotation about the 
axis of symmetry of the spacecraft. Then, we can formulate 
an optimal attitude maneuver problem. 

The performance index is the sum of the magnitudes of the 
initial impulse and the terminal impulse. Equivalently, one 
can minimize the change in the initial angular momentum 
and the change in the terminal angular momentum. Since 
the initial attitude and the terminal time are fixed, Hn and 
Ajv can be considered as functions of no through the discrete 
equations of motion. The optimization problem is equivalent 
to 
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such that C = \\A N - A Nd \\ = 0, 
subject to JgMHl. 

B. Computational approach 

This problem is optimized by the Sequential Quadratic 
Programming (SQP) method using analytical expressions for 
the sensitivity derivatives of the performance index and of 
the constraint equation. 

The variation of the performance index is 
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Since the initial attitude is given and fixed, the perturbation 
of the initial attitude (q is zero. Therefore SH^ — 
and Cjv = ( E>i2^no from J30I . Then, 8 J is given by 
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Since A^ d is fixed and Ajv € S , the variation of the 
constraint can be written as 



6C = -2A T Nd SA N = 2A T Nd R%S(e 3 )C 
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where (jy = 3>i2<5IIo from d30i . Thus, 5C is 

SC = [2Aji?^5(e 3 )$i 2 ] SU . (33) 

Equations (1321 and d33l are analytical expressions for the 
sensitivity derivatives of the performance index and the 
constraint. 

C. Numerical example 

The mass, length and time dimensions are normalized. 
The moment of inertia of the spacecraft is chosen as J = 
diag [3, 3, 2], so that es is the axis of symmetry of the 
spacecraft. 

The desired maneuver is to rotate the axis of symmetry 
from the radial direction to the normal to the orbital plane 
during a quarter orbit. The boundary conditions are given by 

< = diag [1, -1, -1] , A Nd = [0, -1, 0] T . 

We use MATLAB's f mi neon function as an optimization 
tool. The sensitivity derivatives of the performance index 
and the constraint are provided by i32\ and (1331 . The initial 
guess of the initial angular momentum is chosen as iLy = 
J [1, 1, 0] T . The optimized performance index and the vio- 
lation of constraints are J = 6.771, C = 4.80 x 10~ 14 . The 
corresponding angular momenta and the terminal attitude are 

lT 
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-2.915, -2.347, -2.734; 

= [-0.343, 2.686, 2.734] T 

0.633 0.733 0.000 " 
0.000 0.000 -1.000 
-0.773 0.633 0.0000 



so that A N = R%e 3 = [0, 0, -1] T = A^. Fig.|3]shows the 
optimal maneuver of the spacecraft. 

V. Conclusion 

A global model for a rigid spacecraft in a circular orbit 
about a large central body is presented. This model includes 
gravity gradient effects that arise from the non-uniform 
gravity field. 

The sensitivity derivatives for attitude dynamics of a rigid 
body are derived while satisfying the global geometry of the 
problem. Accurate computational approaches for solving a 
nonlinear boundary value problem and the minimal impulse 
optimal control problem for spacecraft attitude maneuvers 
are studied using sensitivity derivatives. 

The attitude dynamics are represented by a rotation matrix 
in the Lie group SO(3), and it is updated by Lie group 
variational integrators that preserve the structure of SO (3) 
as well as other geometric invariants of motion. The sensi- 
tivity derivatives are expressed in terms of the Lie algebra 
so(3). This approach completely avoids the singularities and 
ambiguities associated with Euler angles or quaternions, and 
it leads to a geometrically exact and numerically efficient 
method for rigid body attitude dynamics problems. 

Although the development in this paper includes a gravity 
gradient moment and the rotation of the LVLH frame, the 
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(c) Rotation about the LVLH axes ei and 
Fig. 2. Spacecraft attitude maneuvers 
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Fig. 3. Optimal spacecraft attitude maneuver 



results presented reduce to the case of a free rigid body 
if ujq = 0. That is, the computational approach suggested 
applies directly to attitude maneuvers of the free rigid body. 
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